A resistor model for quantum criticality: quasi one-dimensional organic conductors 
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We investigate the metallic transport in the organic superconductor (TMTSF^PFe under pressure 
within the framework of the spin density wave theory in the proximity of a Peierls quantum critical 
point (QCP). We use a simple transport model of hot and cold regions around the Fermi surface, 
driven by the proximity to the QCP in order to provide a template to fit the experimental results. 
Successful agreement with the data comforts the interpretation of extended criticality around the 
antiferromagnetic QCP in the normal phase of quasi one-dimensional organic superconductors. 

PACS numbers: 74.25.fc, 74.40.Kb, 74.70.Kn 



Anomalous transport properties have been reported 
in the normal phase of the organic superconductors 
(TMTSF) 2 PF 6 and (TMTSF) 2 C10 4 , also called Bech- 
gaard salts [l|, [|| , as a function of temperature and ap- 
plied pressure. These quasi one-dimensional materials 
are remarkable in many respects, the least of all be- 
ing the striking similarity of their temperature-pressure 
(T — P) phase diagram [3| with other exotic supercon- 
ductors, namely the iron-based superconductors 0, Hj]. 
At ambient pressure, (TMTSF) 2 PF 6 undergoes a mag- 
netic phase transition of spin density wave (SDW) char- 
acter, where the Fermi surface is destabilized by the on- 
set of a spin density wave of the itinerant conduction 
electrons [6|-|8(- Under application of high pressure, the 
compound becomes superconducting with a maximum T c 
of 1.2 K at Pqcp ~ 9.4 kbar [l[. This optimal pres- 
sure coincides with the end of the phase of coexistence 
between the antiferromagnetic SDW order and supercon- 
ductivity Q , namely a zero-temperature phase transition 
(QCP), while at larger pressure, superconductivity shows 
a smooth decrease of T c with pressure. The similarity 
with pnictides superconductors is not only seen in the 
(T — P) phase diagram, where one observes a simultane- 
ous occurrence of SDW order and superconductivity, but 
also in the properties of the normal phase Indeed, the 
electrical resistivity p shows a striking linear in tempera- 
ture dependence, p = aT, at the optimal pressure Pqcp 
and above in (TMTSF) 2 PF 6 , suggesting an unusual scat- 
tering mechanism, possibly similar to the one present in 
the pnictides (and cuprates). The same behavior was 
reported in (TMTSF) 2 C10 4 above T c . 

Extensive theoretical investigation has been carried 
out, based on a one-loop renormalization group (RG) for- 
malis m l lOlj especially suitable for one-dimensional sys- 
tems [Tl|. These studies have been reported to repro- 
duce many of the most striking features observed exper- 
imentally [l(J 12 1 . Of special importance is the inter- 



play between itinerant antiferromagnetism (SDW) and 
superconductivity, which controls the magnitude of the 
coupling constants in the vicinity of the QCP. An exten- 



sive regime of criticality has been uncovered, for which 
linear in T resistivity is obtained down to lowest tempera- 
tures (within the one-dimensional framework) on a finite 
range of pressure. Furthermore, the T-linear behavior 
disappears above an upper scale T where the quadratic 
behavior is recovered. Above 9.4 kbar, the system is 
thought to be homogeneous, and although the quantum 
phase transition is reported to be weakly first-order 
it is believed that considerable quantum fluctuations are 
present to cause anomalous regimes, in both transport 
and thermodynamics. 

In this Letter, we build upon the RG insight and sug- 
gest an explicit model based on a Peierls instability in a 
quasi-onc-dimensional system driven to zero-temperature 
by the presence of warping on the Fermi surface. The 
idea of driving a finite temperature Peierls transition to- 
wards a QCP end point through increasing the warping 
is not new and has been used in past studies of chromium 
14 1, a system showing almost perfect nesting be- 
tween Fermi pockets via translation of the SDW wave 
vector. Here, we shall employ a model of the same type, 
where the SDW transition is driven to zero-temperature 
through the increase of the unnested part of the elec- 
tronic dispersion [I3-[l7j]. To the best of our knowledge, 
the theoretical investigation of a Peierls-type QCP driven 
by curvature effects is unprecedented while the family of 
organic superconductors provides a suitable model sys- 
tem to check the theory. 

In the Bechgaard salts, the kinetic electronic spectrum 
is suitably modeled by the orthorhombic dispersion rela- 
tion, 

e(p) = v (\p x \ - p F ) - t b cos(bp y ) - t' b cos(2bp y ) (1) 

The corresponding Fermi surface is shown in Fig. [TJa). 
Generally t' b <C tf, <C vqPf and the t'^-unnesting term 
drives the system into criticality. For t' b = 0, the Fermi 
surface is perfectly nested with the nesting vector Q = 
(2pF, ft/b) and at sufficiently low temperature the system 
is found to be in the SDW phase. For finite yet still 
small t' b , nesting is no longer perfect but still very good in 



the proximity of the inflection points — the hot spots — 
of the Fermi surface. Eventually as t' b is increased further 
(by applying greater external pressure) beyond a critical 
value, nesting will be ineffective and the SDW order is 
destroyed. Nesting of finite parts of the Fermi surface at 
high temperature is now reduced to four nesting points 
at T = as depicted on Fig. QJa) . This picture survives 
down to zero-temperature, implying a quantum phase 
transition at a critical value of the coupling constant t' b . 
Using the parameters of Ref. [l(J, we have %~ 200 K, 
and the critical coupling t' b sa 25.4 K. 

Looking for nontrivial effects such as the observed 
QCP, we need to enrich the model of noninteracting 
fcrmions by a proper model for the two-particle interac- 
tion. Previous works established that there are three rel- 
evant interaction channels in the Bechgaard salts: back- 
ward scattering (with the coupling constant g\ ), forward 
scattering (g 2 ), and Umklapp scattering (g 3 ) [10j. The 
RG studies have shown that the superconducting fluc- 
tuations lead to a drastic decrease of the coupling con- 
stants gi and g 2l due to the interplay of Cooper and 
SDW fluctuations at the nesting points, whereas the 
Umklapp coupling constant 93 remains essentially unaf- 
fected. For the following QCP study, we represent the 
fermion-fermion interaction as mediated by a bosonic 
mode that becomes critical at the QCP. Following the 
insight of Rcf. icj , we retain only the coupling constant 
173 as medium for the electron-paramagnon interaction. 
Our approach allows to analyze both thermodynamic and 
transport properties, but considering that 33 is related to 
the Umklapp processes enables us to focus mainly onto 
the resistivity behavior. We perform calculations on the 
level of the Gaussian approximation, which we may ex- 
pect to yield at least qualitatively the correct physical 
picture. 

In a phcnomcnological low-energy picture, we may as- 
sume that after integrating out all high-energy degrees 
of freedom, the effective interaction is mediated by long- 
wave paramagnon modes. Here, we consider such a 
bosonic mode that transfers a momentum of order Q, 
cf. Fig. [Ha). The coupling of the bosonic modes to the 
electrons generates a self-energy term Hui,q in the boson 
propagator x(i w i Q + Q.)- I n the one-loop approximation, 
ITj,q is given by the polarization bubble [see Fig. H{b)], 
whose relevant nonanalytic part is 

n 9s\q±\ , f (2Mj) 2 l 

with £ q = vqu — b^q 3 ^ + b^q^. In these formulas, v = |v| 
denotes the absolute value of the Fermi velocity at the 
hot spots and gy (q±) is the momentum component par- 
allel (perpendicular) to v. Note that we have expanded 
the spectrum around the hot spots up to the term of 
forth order in q±, which is driven by the t^-unnesting 
term in the dispersion relation ([T]), 64 oc t' b , and is the 
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FIG. 1: (Color online) (a) Brillouin zone with the Fermi 
surface and the nesting vector Q connecting two inflection 
points Pi and P3. (b) Polarization bubble for the param- 
agnon mode, (c) Fermion self-energy. 



lowest order to break perfect nesting around the inflec- 
tion points. The g^-term is active as soon as the crit- 
ical value for t' b is reached. Importantly, the presence 
of the gj_-term prevents the polarization from produc- 
ing mass terms containing logarithms in temperature and 
thus establishes the existence of a QCP — for subcritical 
values of t' b , the effective absence of forth-order g^-terms 
leads to logarithms InT in the bosonic mass so that the 
phase transition towards an SDW state sets in at a tem- 
perature > 0. The remaining analytic part of the bosonic 
spectrum is generically an analytic function of £q. As a 
result, we write the effective propagator for the para- 

magnons as x( iw , Q + q) = [m + a (W w ) 2 + n w , q ] 
with a ~ 1. The bosonic mass \x measures the dis- 
tance to the QCP and, close to it, we should consider the 
limit [i — > 0. The logarithm present in the paramagnon 
propagator is characteristic of a Peicrls phase transition 
and all the anomalous behavior shall ultimately be due 
to this nonanalyticity. 

Using this effective model for the hot-spot electrons 
coupled to critical paramagnons, we are in a position 
to investigate their transport properties. The relevant 
quantity is the retarded electron self-energy T, R (e), which 
within the precision of the one-loop approximation is rep- 
resented by the Feynman diagram in Fig. [He). I 11 the 
standard way for an itinerant electron QCP, the momen- 
tum dependence of the self-energy is negligible compared 
to the energy dependence and the Matsubara self-energy 
is thus given by 

^) = - g3 TE/^ T ^H (3) 
^ J 2nv i(e + u>) - £ q 

with x(iw,£ q ) = J(dgj_/27r) x(iw,Q + q)|<r q=const - Per- 
forming the analytical continuation to real-time frequen- 
cies e, we obtain in the limit e — ¥ for the imaginary 
part of the retarded self-energy in the electronic Green's 
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functions the formula 



Im E R (e) ~ ttT 



In (p F 2 n + e 2 /e 2 F ) 



In 



(4) 



It shows that at criticality, p = 0, the self-energy of 
hot-spot electrons is linear in temperature, ImSq CP (£— > 
0) = 7tT, and independent from the coupling constants. 
As a straightforward consequence, the resistivity of the 
hot-spot electrons in the compound would at arbitrarily 
low temperatures be linear in T as well. Away from the 
QCP, the finite bosonic mass p suppresses for frequen- 
cies |e| < Vy/JI the quantity ImE^e), which for e —> 
tends logarithmically to zero. As a result, since for the 
conductivity essential frequencies e are of order T, a lin- 
ear law for the temperature dependence of the hot-spot 
resistivity appears only above the critical temperature 



T S ~ Vy/Jl . 



(5) 



At the QCP, clearly, T s = 0. Note that, since the 
limit T 3> Tg is essentially equivalent to the limit p — >■ 0, 
the coefficient in front of T in the resistivity does effec- 
tively not depend on the value of the bosonic mass p, i.e. 
on the applied pressure that determines p. 

We turn now towards a dichotomic description of the 
transport properties 18[ of the compound in terms of 
hot-spot and cold-spot regions on the Fermi surface, 
hereby providing a simple model upon which to test the 
experimental data of (TMTSF) 2 PF 6 from Ref. Q. The 
conductivity <j{T) is the sum of contributions from the 
entire Fermi surface. Treating separately the contribu- 
tions from the hot spots [with the Fermi surface volume 
fraction Vh] and those due to the cold regions [volume 
1 — Vh], we write cr(T) as the sum 



a(T) 



1 - v h 



Pa + Phot(T) po + p C oid(T) 



(6) 



In circuit language, see the inset of Fig. O this formula 
corresponds to the parallel arrangement of the resistances 
due to the hot and cold regions of the Fermi surface while 
each of the two resistances is viewed as a series of the 
residual resistance and a specific temperature-dependent 
one. The residual resistivity p is experimentally given 
by the T — > limit and is the result of elastic scat- 
tering processes. Guided by the preceding theoretical 
considerations, we specify in the following the form of 
the temperature-dependent resistivities phot/coid(T) and 
their underlying scattering processes in the hot and cold 
regions. 

For the cold regions, we may for all temperatures as- 
sume the quadratic law p co id(T) = BT 2 accounting for 
the electron-electron scattering processes typical of the 
metallic behavior. At a sufficiently high temperature 
T > To, the notion of cold and hot regions is irrele- 
vant so that we may expect the same law also in the 
hot regions, pi w t(T) = BT 2 . Lowering the temperature, 



TABLE I: Temperature dependencies of the resistivities of hot 
and cold electrons. 



temperature region 


Phot(T) 


Pcold(T) 


T > To 


BT 2 


T s <T < T 


AT 




BT 2 


T <T S 


CT 2 




BT 2 



we encounter at a temperature Tq the crossover into the 
quantum critical regime. Here, (Umklapp) scattering of 
hot conduction electrons through the quantum critical 
paramagnons leads according to the preceding analysis 
to a linear law phot(T) = AT, cf. Eq. (j3J). Below a second 
crossover temperature Ts, Eq. (]5|), the linear resistivity 
is suppressed and one should again expect a Fermi-liquid 
like behavior, /Ohot(T) = CT 2 , though with an effective 
quasi-particle mass heavily renormalized by the interac- 
tion with paramagnons close to criticality. At the QCP, 
Tg = so that the linear law for Phot(T) prevails down to 
zero temperature while at very high pressure, the differ- 
entiation between hot and cold regions is no longer valid 
so that we expect C — > B and the critical window be- 
tween Ts and To to shrink to zero. Table U summarizes 
the temperature laws for the three regimes. 

Applying the dichotomic description (|6]) to the exper- 
imental data on the transport in (TMTSF) 2 PF 6 [!], the 
coefficient B is fixed from the quadratic resistivity law 
po + BT 2 at high temperatures (~ 30 K) while the resid- 
ual resistivity po is given by the T — > intercept of the 
data. Then, we extract the critical regime Tg < T < To 
where the fitting procedure yields according to Eq. ((6]) 
the coefficient A. Last, we use the low-temperature 
(T < Ts) form of Eq. ^ to determine the coefficient C. 
Within the philosophy of the resistor model, A and C 
should be kept constant function of temperature 

inside the temperature regime they appear in. This en- 
sures that the regime is properly defined according to 
Eq. © and Table fl] Theoretically, we may expect log- 
arithmic corrections, see Eq. Q, but when compar- 
ing with experiments, these are fairly approximated by 
constants. Finally, we determine the crossover temper- 
atures Tg and To as the intersections of the asymptotic 
forms for each regime. 

Figure [2] shows the pressure dependence of B, C, and 
Vh as a result of the fitting analysis of <r(T) between 0.15 
and 34 K at seven different pressures according to the 
procedure discussed above. In line with the theoretical 
prediction, we may assume that A is independent from 
pressure. Treated as a free fitting parameter, it would 
mildly jitter around A = 0.38 pfl cm/K, confirming the 
assumption of A being constant. The coefficient B is re- 
lated to the effective mass in the cold regions, B ~ mf,. 
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FIG. 2: (Color online) Pressure dependence of the model pa- 
rameters obtained from the fit: B (diamond, in fifl cm/K 2 ), C 
(square, in fiQ cm/K 2 ), and vh (circle) compared to the tem- 
perature of the superconducting transition T c (triangle, in K) 
from Ref. Q. Inset: "equivalent circuit" for the dichotomic 
conductivity model ©. 
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FIG. 3: (Color online) (pressure— temperature) phase dia- 
gram of (TMTSF^PFe displaying the crossover tempera- 
tures T$ and To from this analysis as well as the long-range 
order phases. The lines are guides to the eye. 



Its slight decrease under pressure in Fig. [2] can be as- 
cribed to the increase of the intermolecular in-chain over- 
lap, possibly enhanced by correlations. The data are also 
in agreement with the pressure dependence of the spin 
susceptibility measured by NMR experiments [l9[. The 
coefficient C describing the increase of the effective elec- 
tron mass at hot spots is roughly constant but its size is 
ten times larger than the order of magnitude of B. The 
fact that such an enhancement does not fade away under 
pressure indicates that even under 21.8 kbar the scat- 
tering off antiferromagnetic spin fluctuations is still very 
strong. Under a pressure of 11.8 kbar, corresponding to 
the point closest to the QCP, no quadratic law could be 
observed down to the lowest temperature after supercon- 
ductivity had been removed by the application of a small 
magnetic field along c*. This suggests that the quantum 
critical regime of linear resistivity extends down to tem- 
peratures very close to zero at this point. The hot-spot 
volume Vh, which close to the QCP is « 0.97, is de- 
creasing under pressure (v^ sa 0.30 at 21.8 kbar). This 
is in accordance with the intuitive physical picture that 
the distance from the QCP enhances unnesting and thus 
reduces the effective size of the hot spots. Its value re- 
markably follows T c , in agreement with earlier findings 
I- 

Both crossovers To and Tg are plotted versus pressure 
in Fig. |3] The behavior of Ts is strongly suppressed in 
the vicinity of the QCP, in fair agreement with Eq. (0. 
It is to be noted that while the hot-spot contribution 
to conductivity is dominant at the pressure of 11.8 kbar 
close to the QCP, the presence of the cold regions is cru- 
cial to explain the pressure decrease of the resistivity at 
a fixed temperature. Indeed, as discussed above, increas- 
ing the distance to the QCP induces a decrease of Vh, 
thus favoring the conduction through the cold regions at 
larger pressures. In the language of the equivalent cir- 



cuit (see inset of Fig. [2]), the less-resistive cold regions 
short-circuit the larger and for Tg < T < T linear in T 
hot-spot resistance. 

In conclusion, we present a theory of a QCP associated 
to the Peierls-type singularity. This theory is nontrivial 
as the role of the curvature is preponderant in stabilizing 
the logarithmic divergences and yields strong influence 
on the form of the crossovers. The organic Bechgaard 
salts constitute an almost perfect model system with the 
simplicity of their band structure allowing to test the cur- 
vature effects. Within a hot-spot/cold-spot dichotomic 
conductivity model for the itinerant electron QCP, we 
confront the critical theory with the experimental data 
obtained in transport measurements for (TMTSF^PFg, 
showing a good agreement. At the hot spots, the physics 
of the Bechgaard salts shows strong similarity with the 
physics in heavy-fermion systems. 
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